ON THE QUALITATIVE BEHAVIOUR OF INCOMPRESSIBLE 
TWO-PHASE FLOWS WITH PHASE TRANSITIONS: 
THE CASE OF EQUAL DENSITIES 

JAN PRUSS, GIEPJ SIMONETT, AND RICO ZACHER 



Abstract. The study of the basic model for incompressible two-phase flows 
with phase transitions in the case of equal densities, initiated in the paper 
Priiss, Shibata, Shimizu, and Simonett I16| , is continued here with a stability 
analysis of equilibria and results on asymptotic behaviour of global solutions. 
The results parallel those for the thermodynamically consistent Stefan prob- 
lem with surface tension obtained in Priiss, Simonett, and Zacher |19| . 
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1. Introduction 

In this paper we study a sharp interface model for two-phase flows with surface 
tension undergoing phase transitions. The model is based on conservation of mass, 
momentum and energy, and hence is physically exact. It further employs the 
standard constitutive law of Newton for the stress tensor, Fourier's law for heat 
conduction, and it is thermodynamically consistent. 

Suppose that two fluids, fluidi and fluid2, occupy the regions Qi(t) and ^(i), 
respectively, with U SI2 (t) = ^- Let T(t) = dfli(t) be a sharp interface that 

separates the fluids. Across the interface T(t) certain physical parameters, such as 
the density, viscosity, heat capacity and the heat conductivity, experience jumps. 
We assume that the interface is ideal in the sense that it is immaterial, which 
means that it has no capacity for mass or energy except surface tension. 

In more detail, let C K™ be a bounded domain of class C 3_ with n > 2. We 
further assume that T(t) n dfl = 0, which implies that no boundary contact can 
occur. In the following we let 

• Ui denote the velocity field in fli(t), 
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• TTi the pressure field in fij(f), 

• Ti the stress tensor in fii(i), 

• Di — (Vtti + [Vw,i] T )/2 the rate of strain tensor in flt(i), 

• Qi the (absolute) temperature field in f2j(t), 

• the outer normal of Q,\(t), 

• up the velocity field of T(t), 

• Vr = ur • the normal velocity of T(t), 

• i?r = H(T(t)) = — divr^r the sum of the principal curvatures of T(t), and 

• = V2 — v i the jump of a quantity u across r(t). 

Here the sign of the curvature Hp is negative at a point x £ T if f2i n B r (x) is 
convex, for some sufficiently small r > 0. Thus if fii is a ball, i.e. T = Sr(xq), 
then H T = -(n-l)/R. 

Several quantities are derived from the specific free energies tpi(9) as follows: 

• Ci(9) = ipi(9) + 0i]i(9) is the internal energy in phase i. 

• T)i(9) = —ip'i{6) is the entropy, 

• Ki(9) = e'^9) = -6tp"(9) > is the heat capacity, 

• 1(9) = 9{i)'(9)j = -9{r](6)] is the latent heat. 

Further di(9) > denotes the coefficient of heat conduction in Fourier's law, 
fJ,i(9) > the viscosity in Newton's law, p := pi = p<i = 1 the constant density, 
and a > the constant coefficient of surface tension. 

In the sequel we drop the index i, as there is no danger of confusion; we just 
keep in mind that the physical quantities depend on the phases. 

By the Incompressible two-phase flow with phase transition we mean the fol- 
lowing free boundary problem: find a family of closed compact hypersurfaces 
{r(i)} t >o contained in f2 and appropriately smooth functions u : R + x f2 — > R™, 
and 7r, 9 : M+ xii->l such that 





dtu + u ■ Vw — div T 


= 


in 


n\r(t) 


T -- 


= p(9)(Vu + [Vw] T ) - n I, divu 


= 


in 


n\r(*) 


n(9)(d t 9- 


f u ■ V9) - div(d(9)V9) — T : Vit 


= 


in 


n\r(t) 




u = d v 9 


= 


on 


on 


< 


m = m 


= 


on 


r(t) 




[Tur] + o-Ht-vt 


= 


on 


r(t) 




\i){9)\ + aH T 


= 


on 


r(t) 




-ld{9)d v 9}+l{9){V v -u-u r ) 


= 


on 


r(t) 


r(o) 


= T , u(0,x) = uo(x), 9(0, x) 


= 9 (x) 


in 


Q. 



This model has been recently proposed by Anderson et al. p], see also the mono- 
graphs by Ishii |9 j and Ishii and Takashi [TU] , and the derivation in Section 2 of the 
recent paper [16] . It has been shown in |16j that the model is thermodynamically 
consistent in the sense that in the absence of exterior forces and external heat 
sources, the total energy is preserved and the total entropy is nondecreasing. It 
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is in some sense the simplest sharp interface model for incompressible Newtonian 
two-phase flows taking into account phase transitions driven by temperature. 

There is a large literature on isothermal incompressible Newtonian two-phase 
flows without phase transitions, and also on the two-phase Stefan problem with 
surface tension modeling temperature driven phase transitions. On the other hand, 
mathematical work on two-phase flow problems including phase transitions are 
rare. In this direction, we only know the papers by Hoffmann and Starovoitov 
|H] dealing with a simplified two-phase flow model, and Kusaka and Tani [T31 [2] 
which is two-phase for temperature but only one phase is moving. The papers of 
DiBcnedetto and Friedman [2] and DiBenedetto and 0'Leary[3] deal with weak 
solutions of conduction-convection problems with phase change. However, none of 
these papers deals with models which are consistent with thermodynamics. 

It is the purpose of this paper to present a qualitative analysis of problem 
(|l.lj) in the framework of i p -theory. We discuss the induced local semiflow and 
study the stability properties of the equilibria. These are the same as those for the 
thcrmodynamically consistent two-phase Stefan problem with surface tension, and 
even more, also their stability properties turn out to be the same. This heavily 
depends on the fact that the densities of the two phases are assumed to be equal; 
in this case the problem is temperature dominated. 

In a forthcoming paper we will consider the case where the densities are not 
equal; then the solution behavior is different, as the interfacial mass flux has a 
direct impact on the velocity field of the fluid, inducing so-called Stefan currents. 
The velocity field is no longer continuous across the interface which leads to dif- 
ferent analytic properties of the model. We call this case velocity dominated. 

It has been shown in [TB] that the total energy 

E := E(u, 9, r) := i f \u\ldx+[ e{9) dx + cr\T\, (1.2) 
2 Jn\r Jn\r 

is preserved along smooth solutions, while the total entropy 

$(0,T) = ( n{9)dx (1.3) 

is strictly increasing along nonconstant smooth solutions. By similar arguments as 
in [19] , it can further be shown that the equilibria of (jl.ll) are precisely the critical 
points of the entropy functional with prescribed energy, and that a necessary 
condition for such a point e* = (u*, r») to be a local maximum of the entropy 
functional with prescribed energy is that T* is connected and that the stability 
condition (S), see Theorem 13.11 below . is satisfied. 

The plan for this paper - which builds on [TB] and [TO] - is as follows. Our 
approach is based on the so-called direct mapping method where the problem with 
moving interface is transformed to a problem with fixed domain, resulting in a 
quasilinear parabolic evolution problem with a dynamic boundary condition on 
a domain with fixed interface. The main result on well-posedness of the trans- 
formed problem is taken from [16) and is stated in Section 2. The linear stability 
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properties of the equilibria are derived in Section 3. It turns out that generically 
the equilibria are normally hyperbolic. They are always unstable if the disperse 
phase f^i is not connected. If both phases are connected we find the same stability 
condition, condition (S) in Theorem 13. II below, as in Priiss, Simonctt and Zacher 
[19], see also Priiss and Simonett [17] • As the equilibria are normally hyperbolic 
we may use a variant of the generalized principle of linearized stability, see Priiss, 
Simonett and Zacher |18j . to prove nonlinear stability or instability. Combining 
this method with the Lyapunov functional we are able to show that a solution 
which does not develop singularities exists globally and its orbit is relatively com- 
pact in the state manifold. If such a solution contains a stable equilibrium in its 
limit set, then it is shown that it converges to this equilibrium. 

2. The Local Semiflow 

(i) Local Existence 

The basic result for local well-posedness of problem (jl.lj) in an L p -setting, stated 
in [16] Theorem 5.1], is the following. Here Pr = I — vy®vt denotes the orthogonal 
projection onto the tangent space of T. 

Theorem 2.1. Let p > n + 2, a > 0. Suppose tpi € C 3 (0, oo), fi u d l € C 2 (0, oo) 
such that 

«i(s) = -sip"(s) > 0, (j,i(s) > 0, di(s)>0, s e (0,oo), i = 1,2. 
Assume the regularity conditions 

(u ,8 ) e [w;^(fi\r )nc(n)r +1 , r e w^, 

where fl C M™ is a bounded domain with boundary dfl £ C 3 ~ , 
the compatibility conditions 

divuo = in Q \ To, uq — <9„#o — on dQ, 

ProIM^oXVuo + [V Uo ] T )l = on T , 

ty(uo)] + <JH F0 = on T , [d^oJS^Uo] e W^ 6 'p{T ), 

and the well-posedness condition 

6q > on f2, 1(uq) ^0 on Tq. 

Then there exists a unique L p -solution of problem (| 1 . 1|) on some possibly small 
but nontrivial time interval J = [0, r] . 

(ii) The Local Semiflow 

We follow here the approach introduced in Kohne, Priiss and Wilke [TT] for the 
isothermal incompressible two-phase Navier-Stokes problem without phase transi- 
tions and in Priiss, Simonett and Zacher [19] for the Stefan problem with surface 
tension. 
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Recall that the closed C 2 -hypersurfaces contained in fi form a C 2 -manifold, 
which we denote by MH 2 {^1). The charts are the parameterizations over a given 
hypersurface E, and the tangent space consists of the normal vector fields on E. 
We define a metric on AiH 2 ^) by means of 

where dn denotes the Hausdorff metric on the compact subsets of K™ and A/" 2 E = 
{p, vz(p), Ve&"e(p)) : peS} the second order bundle of E € A4% 2 (f2). This way 
■M"H 2 (S7) becomes a Banach manifold of class C 2 . 

As an ambient space for the state manifold SM of problem (II. ip we consider 
the product space C(Cl) n+1 x VWH 2 (f2), due to continuity of velocity, temperature 
and curvature. 

We then define the state manifold SM as follows. 

SM := {(u,e,T) G C(0)™ +1 x MU 2 : (u,d) G w£- 2 / p (fi \ r) n+1 , r G W^ 3/p , 

div it = in D,, 9 > in CI, u = d v 9 = on 90, 

PrK^^r] = 0, [ip(6)l +aH r = on T, (2.1) 

l(9)^0onT, \dd u 9\ G ^ 2 - 6 ^(r)}, 

Charts for these manifolds are obtained by the charts induced by VWH 2 (f2), fol- 
lowed by a Hanzawa transformation. 

Applying Theorem 12.11 and re-parameterizing the interface repeatedly, we see 
that (jl.ip yields a local semiflow on SM. 

Theorem 2.2. Letp > n+2, a > 0, and suppose ipi G C 3 (0,oo), <ij G C 2 (0,oo) 
such that 

Ki(s) = -sip"(s) > Q, m(s) > 0, di(s)>0, s G (0,oo), i = 1,2. 

T/ien problem generates a local semiflow on the state manifold SM. Each 

solution (it, 9, r) exists on a maximal time interval [0, £»), where i* = f*(wo, #0j To). 

Note that the pressure does not occur explicitly as a variable in the local semi- 
flow, as the latter is only formulated in terms of the temperature 9, the velocity 
field u, and the free boundary T. The pressure n is determined for each time t 
from (u, 9, T) by means of the weak transmission problem 

(W|V^) L2(fi) = (2div( M (0)£>) - u ■ v«|v^) i8(n)J ^ G flj(n), 

[tt] = o-iJr + 2\n(ff)Dvr • it] on I\ 

Concerning such transmission problems we refer to [111 Scetion 8]. 

3. Linear Stability of Equilibria 

1. As shown in [Tj)J Section 3], the equilibria (u*, 7r „,,#», T*) of (jl.ip consist of 
zero velocities u*, constant pressures 7r* in the phases, constant temperatures 9*, 
and fii is a ball Cl\ — B&,(a:*) C 51 in case f2i is connected, and a union of 
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nonintersecting balls of equal radii otherwise. We assume here that the balls do 
not touch the outer boundary <9f2, to avoid the contact angle problem, and we 
also assume that the balls do not touch each other. We are not able to handle 
the latter case as the interface T* = dfli will then not be a C 2 -manifold. We call 
such equilibria non-degenerate. The temperature 9* and the pressure jump [Vr*] 
are related to i?* via the curvature Hr t through the relation 

(n - l)er 



-<rH r , 



-M9* 



(3.1) 



In the sequel we only consider non-degenerate equilibria and denote the set of such 
equilibria by £ , i.e., 

£ = j(o,0„r„):r„= {Jrt i1 = Sfl,(^)j, 

with [7r*], 9 ^ and i?* determined by (|3.1j) . According to (11.21) the total energy at 
an equilibrium (0,#*,r*) is then given by 



<f{9* 



£(o,0*,r*) 



e(6**) dx + fj|r„ 



n\r» 



(3.2) 



By employing the Hanzawa transformation, see [161 Section 2] , one shows that 
the fully linearized problem at an equilibrium is given by 



dtu — /j,*Au + V7r 


fu 


in 




divu 


= Id 


in 




u 


= 


on 




M 


= 


on 


r* 


-[Tv\ + a(A«h)v 


= 9u 


on 


r* 


K*3ti? — d*A$ 


= .fe 


in 






= 


on 


on 


¥1 


= 


on 


r* 


l*$ — oA*h 


= 9e 


on 


r* 


(h/9*)(d t h-u-v)-ldM] 


= 9h 


on 


r* 


u(0) = uo, 0(O)=i?o, HO) 


= V 







(3.3) 



with 

/z, = M(0*), k* = k(6»*), I* = 1(9*), A* = -H'(0) = -(n - 1)/-R» 



Ar„ (3.4) 
— 6 '*) / '6 '* is 



where Ar» denotes the Laplace-Beltrami operator of r», and i9 = 
the relative temperature. 

It follows from the maximal regularity results in 16 that the operator L defined 
by the left hand side of (|3.3|) is an isomorphism from E( J) into R(J) C F(J) x 
7E, where K(J) is determined by the natural compatibility conditions. Here the 
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function spaces E(J), 7E(J) and F(J), with J = [0, a] an interval, are defined as 
follows: 

E(J) := {(u,7r,g,0,/i) G Ei (J) x E 2 (J) x E 3 (J) x E 4 (J) x E 5 (J) : g= [tt]}, 
where 

Ei (J) :={ U ei7 p 1 (J;L p (r!)) n nL p (J;i/2( r! \ r ^)n :u = 0onffi, [u] = 0}, 
E 2 (J) :=L p (J;ij p 1 (r!\r,)), 

E 3 (J) := W p 1 / 2 -V2p(j ; Lp (T,)) n L p (J; W p 1_1 / p (r*)), 

E 4 (J) :={0eJ^(J;Lp(fi))nL p (J;fl2(fi\r.)) : 9,0 = on 50, [0] = 0}, 
E 5 (J) := W p / 2_1/,2p (J; L P (T,)) n W^ 1 " 1 ^./; ff p 2 (r*)) n L P (J; W^ 4 " 1 ^ (T,)). 
The time- trace space 7E( J) is given by 

7 E(j) = {(«o,i?o,M e (T^ 2 - 2 /f(fi\r,)nc(f7))" +1 x w^(T,)}, 

while the space of data is F(J) := {{f u , fd,9u, fe,96,9h) G ]lj=i F j(-0}> where 
Fi(J) :=i p (J;i p (0))", 

F 2 (J) := ^(J; ff- 1 ^)) n L p (J; H^Q \ r*)), 

F 3 (J) := W^-^i J; L p (r,)) n n L p ( J; W; 1 " 1 ^ (r,)) n , 
F 4 (J) :=L p (J;L p (Q)), 

F 5 (J) := ^^/^(J; L P (I\.)) n L p (J; W 2 " 1 ^ (r.)), 
F 6 (J) := W p 1 / 2 -V2 P(J; Lp(I%)) n ip(J; Wp i-VP (r ,)). 

If the time derivatives d t are replaced by d t + u, with u > sufficiently large, 
then this result is also true for J = M + . 

2. We introduce a functional analytic setting as follows. Set 
X = L p , a (Q) n x L p (fi) x ^ 2 - 2 /f(r»), 
where the subscript cr means solenoidal, and define the operator L by 

L(u,6,h) = ( - y*Au + V7r,-(d»/K»)Ai?,-u- 1/ - (fl./Z.)^*^]). 
To define the domain -D(L) of i, we set 

Xi = {(u,0,/i) g (ff p 2 (ft\r,)nc(fi)) n+1 x wf-V'ir.) ■. 

divu = in u = d u -d = on <9fi}, 

and 

D(L) = {(u,0,/i) G Xi : [P»/i»£>i/] =0, = on r», 

[d„a,0] g w p 2 - 2 /f(r,)}, 



s 
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where P* = Pr„ denotes the orthogonal projection onto the tangent space of T*. 
Here tt is determined as the solution of the weak transmission problem 

(vtt|v0) 2 = (m*au|V0) 2 , 4> e ^>(n), 

[tt] = -crA»/i + 2[/i*£)i/-z/]. 

We refer to [UJ Section 8] for a detailed analysis of such transmission problems. 
The linearized problem can be rewritten as an abstract evolution problem in Xq, 



z + Lz = f, t>0, z(0) = zq, 



(3.5) 



where z = (u, &, h), f = (f u ,fe,9h), z o = ( u o, $o, ho), provided f d = g u = 9e = 0. 
As the terms u ■ v and aA^h are of lower order we may deduce maximal ir- 
regularity of (|3.5|) from that of the decoupled system (cf. (TTJ Section 6] and [EJ 
Theorem 4.3]) by a perturbation argument. We can then conclude that — L gen- 
erates an analytic Co-semigroup in Xo; cf. Priiss [151 Proposition 1.1]. 

3. The eigenvalue problem. Since the embedding X\ Xo is compact, the 
semigroup e~ Lt as well as the resolvent (A + L)^ 1 of — L are compact. Therefore, 
the spectrum a(L) of L consists only of countably many eigenvalues of finite al- 
gebraic multiplicity and is independent of p € (l,oo). Therefore it is enough to 
consider the case p = 2. In the following, we will use the notation 



(u\v)n := (u\v)l 2 (Q) ■= / uvdx, u,v G L 2 {yi) 



(g\h)r, '•= (g\h) L2 (r,)'-= J ghds, g,heL 2 (T«), 

for the L2 inner product in Q and r*, respectively. Moreover, we set \v\q — (v\v 
1/2 

and \g\r, — (<?|<?)r ■ The eigenvalue problem for — L reads as follows: 



1/2 



Xu — /l,Au + V7T 


= 


in 


n\r, 


divw 


= 


in 




u 


= 


on 


dfl 


H 


= 


on 


1 


-\Tv\+<T{A.h)v 


= 


on 


r*, 



(3.6) 





- d*Ad 


= 


iii 


n\r* 






= 


on 


on 




m 


= 


on 


r* 




— aA*h 


= 


on 


r» 






= 


on 


r*. 



(3.7) 
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We are now ready to formulate the main result of this section. 

Theorem 3.1. Let L denote the linearization at (0, #*,r*) G £ as defined above. 
Suppose I* =/= 0. Then —L generates a compact analytic Co-semigroup in Xq which 
has maximal L p -regularity. The spectrum of L consists only of eigenvalues of finite 
algebraic multiplicity. Moreover, the following assertions are valid. 

(i) If r* is connected and the stability condition 

(B) .:-*J:_£££i>- «™ S 
R% 9*{k*\1)q 

holds, then all eigenvalues X ^= of —L have negative real part. 

(ii) The stability condition (S) is equivalent to (//(#*) < 0, where the function 
ip is defined in (|3.2[) . 

(iii) // r* = U™ =1 T^ and s > 0, then —L has precisely m positive eigenvalues, 
and precisely (m — 1) positive eigenvalues if s < 0. 

(iv) X — is an eigenvalue of L with geometric multiplicity (mn + 1). It is 
semi-simple if s =/= 0. 

(v) Let e* = (0, 0*, T*) G £ be an equilibrium. Then in a neighborhood of e* the 
set of equilibria £ forms a (mn + 1)- dimensional C 1 -manifold. Moreover, 
the kernel N(L) of L is isomorphic to the tangent space T er £ of £ at e*. 

Consequently, (0,9*, T*) G £ is normally stable if and only if s < and T* is 
connected, and normally hyperbolic if and only if s > 0, or T* is disconnected and 

S y£ 0. 



Proof, (i) Suppose A with Re A > is an eigenvalue of — L with eigenfunction 
(it,»9, h). Taking the inner product of the eigenvalue problem (I3.6[) with u and 
integrating over 51 we get 

= A|u|n - (div T\u)n = X\u\q + / T:Vudx+ / (T 2 vu 2 - T x vui) ds 

= X\u\ 2 u + 2\ f /J 2 D\l + (pV]| u ) r . 

= X\u\l + 2\nl /2 D\ 2 n + aX(A*h\h) r , + a(A if h\j) v , , 

since \u\ — 0, \Tv\ — a(A*h)v and u ■ v = Xh + j with (U/9*)j = — [d* On 
the other hand, the inner product of the equation for i? with d and an integration 
by parts leads to 

= X\nJ 2 d\l + |di /2 W|2 + ([d*d„tfp) r , 
= X\kI /2 ^\ 2 2 + \dl /2 V$\l - *V\A,h) T J6„ 
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where we employed the relations (Z*/0*)j = — [d*c)„'#] and = aA*h. Adding 
these identities and taking real parts yields the important relation 

= ReAld^ + 2\[il /2 D\l + crRe A (AMh) r 

(3-8) 

in • 



0* (ReA| K i /2 i?| n + \dl /2 Vtf\ 2 



On the other hand, if ft :— ImA 7^ 0, then taking imaginary parts separately we 
get with a = u{A if h\j)Y, 

= Plufa - a/3(A*h\h)r, +lma, 
= PO^kJ 2 ®^ +lma. 

Hence 

a(A*h\h) T , = \u\ 2 n -9*\Kl /2 d\ 2 n . 
Inserting this identity into (|3.8p leads to 



= 2Re \\u\l + 2\[il /2 D\ 2 n + 8*\dl /2 V0\ 2 Q , 

which by (|3.6p - (|3.7[) (and Korn's inequality in case Re A = 0) shows that if A is an 
eigenvalue of — L with Re A > then A is real. 

Supposing that A > is an eigenvalue, we decompose i9 = i?o + i9, h = Iiq + h 
and j = jo + j , where 

£ = (K.|0W(K.|l)n, h = (h\l) r J\T m \, j = (j|l)r»/|r*| 
are weighted means. Then 

\ K l /2 tf\ 2 n = \ K l /2 l) \ 2 n +( K *\l)n$ 2 , \h\l t = \h \ 2 r ,+ \r*\h 2 , blr, =|jolr,+ |r*|j 2 . 
Therefore ()3.8j) becomes 

= \\u\l + 2\ f il /2 D\l + *\(A+ho\ho) r , 

+ 9, (\\ K l /2 \l + \dl /2 VMl) + \6*(K*\l) n ti 2 - \a^-\T*\h 2 . ^ 
We further have 

X J hds = J (u ■ v — j) ds = — J j ds 
hence Ah = — j. Also, the identity 

(l*/6*) [ jds = - [ {dMjds= [ d*A$dx = \ [ K*-&dx 
Jr, Jr, Jn Jn 

implies (Z*/0*)|r*|j — A(K*|l)n^. Thus (|3.9p becomes 
= A|u| n + 2\nl /2 D\l + a\{A*h \h ) rt 



. (Almoin + \dl /2 VMl)+M^*\{ ] 



l*|r.| <n-l)y h2 (3-10) 



».(k*|i)o Rl 
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As A„ is positive semidefinite on functions with mean zero if T* is connected, in 
this case —L has no positive eigenvalues if the stability condition 

1 2 JT*\ (n-lV 



> 



(3.11) 



0*(>C*|1)2 Rl 

is satisfied. This is the same condition we found for the thermodynamically con- 
sistent Stefan problem with surface tension; see [19] and [17] . 

(ii) The assertion follows immediately from the results in 16, Section 3]. 

(iii) On the other hand, if the stability condition does not hold or if T* is discon- 
nected, then there is always a positive eigenvalue. To prove this we proceed as 
follows. Solve the Stokes problem 



(3.12) 



' Alt — /x* Aw, + V-7T 


= 


in 


n\r. 


div u 


= 


in 




u 


= 


on 




M 


= 


on 


r,„ 




= 9V 


on 


r„. 



and define the Neumann-to-Dirichlet operator 



for the Stokes problem by 



in 

on 
on 



n\r, 

r* 
r 4 



(3.13) 



N^g := u ■ v. Similarly, solve the heat problem 

n*Xd - d*Ati = ' 
d v $ = 

-ld*d v ti\=g 

to obtain d = g, where denotes the Neumann-to Dirichlet operator for the 
heat problem. In the following, we use the same notation for i? and its restriction 

to r*. 

Suppose that A > is an eigenvalue with eigenfunction (u, i9, h). Choosing 
g = —aA^h in (|3.12p we obtain u ■ v = —N^aA^h. Next we solve the heat problem 
(|3.13[) with g = {U/9*){u ■ v — Xh), yielding 

= -(l*/9*)N?(N x s crA*h + Xh). 

This implies with the linearized Gibbs-Thomson law = oA*h the relationship 

-(l 2 J9«)N?(N^o-A*h + Xh) = oA*h, 

hence 

Xh + [JVf + ((£/0»)iVf y^o-A^h = 0. 

Setting 

Tx:= [Ng + W/ejN?)- 1 ]- 1 

we arrive at the equation 

B x h := \T x h + aA*h = 0. (3.14) 
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A > is an eigenvalue of — L if and only if (|3.14[) admits a nontrivial solution. We 
consider this problem in Ii2(T*). Then A* is selfadjoint and 



<r(A*h\h) r . > 



On the other hand, we will see below that and are selfadjoint and positive 
semidefinite on I/2(r*) and hence T\ is selfadjoint and positive semidefinite as 
well. Moreover, since A* has compact resolvent, the operator B\ has compact 
resolvent as well, for each A > 0. Therefore the spectrum of B\ consists only of 
eigenvalues which, in addition, are real. We intend to prove that in case either T* 
is disconnected or the stability condition does not hold, B\ has as an eigenvalue, 
for some Ao > 0. 

To proceed we need properties of the relevant Neumann-to-Dirichlet operators. 

Proposition 3.2. The Neumann-to-Dirichlet operator for the Stokes problem 
(|3.12l) has the following properties in X 2 (r >t ). 

(i) If u denotes the solution of (|3.12[) . then 

(N^g\g)r, = X\u\ 2 n + 2 f ^\D\\ dx, g € L 2 (T„), A > 0. 
Jn 

(ii) For each a G (0, 1/2) there is a constant C > such that 

(N?g\g) r * > ^-^VM., 9 e L a (r,), A > 0. 

In particular, 

|iVf| B( z, 2( r.)) < A >°- 

(iii) Let denote the components of I\ and let be the function which is one 
on r*f } zero elsewhere. Then (N^ g\ek)r, — for each k; in particular N^g has 
mean value zero for each g g L2<T*)- Moreover, with e = J^k e fe we nave ^x e = ® 
and (N$g\e) L2( r t) = for all g E L 2 (i\). 

Proof. The first assertion follows from the divergence theorem. The second as- 
sertion is a consequence of trace and interpolation theory, combined with Korn's 
inequality. The last assertion is implied with div u — by the divergence theo- 
rem. □ 

Proposition 3.3. The Neumann-to-Dirichlet operator for the diffusion prob- 
lem (|3. 13[) has the following properties in Z/ 2 (r„). 

(i) If d denotes the solution of (|3.13[) . then 

{N^g\g)v. = M>fe*\h + \y/d*W\h, g e L 2 (T„), A > 0. 

(ii) For each a € (0, 1/2) and Ao > there is a constant C > such that 

(N?g\g) T . > ^r\N?g\* t , g e L 2 (r„), A > A . 
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In particular, is injective, and 

C 

\N\ |B(i 2 (r»)) < j^i A > A . 
(iii) On 1/2, o(r*) = {<? G L2(r*) : (.a|e)r, = 0}, we even have 

(N?g\g)r, > ^-^Vf 5 e i 2 ,o(r.), A > 0, 

and 

C 

l-Wf lg(^2, (r.)) < ( 1 + A w A>0 - 

In particular, for A — 0, p. 131) is solvable if and only if (g|e)r, = 0, and £/ien t/ie 
solution is unique up to a constant. 

Proof. The first assertion follows from the divergence theorem. The second and 
third assertions are consequences of trace and interpolation theory, combined with 
Poincare's inequality. The last assertion is a standard statement in the theory of 
elliptic transmission problems. □ 

(a) Consider v x := T x e, or equivalently e = N^v\ + (c* A^) -1 i;a, where we used 
the abbreviation c* = l\j and where e is the characteristic function on T*. Here 
I\ can be either connected or disconnected. Denoting the orthogonal projection 
from L2(r*) to £2,o(r*) by Qq, the equation for v\ is equivalent to 

v\ + c.N?Q Q Ngvx = c*iVf e, 

due to Proposition 13.21 Multiplying this identity in ^(r*) by N^v x we obtain 
with Propositions 13.21 and 13.31 

c(A)|JVftfc|£. < (v x + c*N? N^v x \N^v x ) rt = (c*iVf ejiVf v x )r, 

= c,{e\N? Q N^v x )r. < C(A)|iV A 5 i; A |r, , 

where c(A) and C(A) are bounded near A = 0, showing that N^v x is bounded near 
A = 0. This implies 

lim XT x e = lim Xv x = c* lim AiV?e, 

A^O A^O A^O 

provided the latter limit exists. 

To compute this limit, we proceed as follows. First we solve the problem 

— d*A$ = — K*ao in f2 \ T* 
d v ® = on dVL 

m = o on r, (3 - 15) 

-|<i*9 I/ '!?] = e on T*. 

where ciq = \r*\/(n*\l)a, which is solvable since the necessary compatibility con- 
dition holds. We denote the solution by i?o and normalize it by (k*|i?o)o = 0. 
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Then i9a = e — -do — ao/ A solves the problem 

k*Ai? — c?* Ai? = — k^Ai^o in \ 
d v d = on dn 

[ff\ = o, on r„ (3 - 16) 
-[<*.$,0] = o on r„. 

By the normalization («;*|i?o)fi = we see that $ A is bounded in iJf (f2 \ r*) as 
A — > 0. Hence we have 

lim ATVf e = lim [Ai? A + Ai? + a ] = a = |r,|/(«,|l) n . 

A— >0 A->0 



This then implies 



lim( J B A e|e)r t = c*-^— 



if the stability condition does not hold. 

(b) Next suppose that T* is disconnected, i.e. I\ = U^LjTj, and set g = 

X^fc a k^k 7^ with J2k a k = 0. Hence Qog = g. Then for v\ := T A e we have 
as in (a) boundedness of N^v\ and then 

lim XT\q — lim Xv\ = c* lim XNp Qnq = 0, 

A^O A-s-0 A^O 

since N?Q$ is bounded as A — > 0. This implies 



\im(B x g\g) L2{r ,) 



a(n-l) 



< 0. 



(c) Next we consider the behavior of (B\g\g) L2 ( Tt - ) as A — > oo. With c* = 
as above we first have 

T x = (I + c*iVf A^-^JVf = c*JVf - c*iVf JVf (I + c»JVf JVf J-^JVf , 

hence by Propositions 13. 3[ 13.21 for A > Ao, with Ao sufficiently large, 

(lAS|s)r. = c.(JVf s|s)r. - ^(7V A S (7 + c*ATf N^N^g^ g) r . 

l^fU 2( r,) 



> c* 



> c* 



> c* 



(JVfsl«/)r. - c 
«.9l.9)r t - 



\N?g\l 



l-c|JVf| ia( r.)I^U a (r.) 
CVWk(r.) 



l-c|^U.(r.)I^U a (r.) 



(^ffflff)r. 



(JVffll5)r. - 2(^.9l5)r,J = c (ATf fl | 5 ) r . 



Therefore, it is sufficient to bound (N^ g\g)r, from below as A — > oo. 
For this purpose we introduce the projections P and Q by 

m 

Pg = c m ^2(g\e k ) rt e k , Q = I - P, 

k=l 
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where c rn = m/|r*| in case I\ has m components. Then with g k — (ff|efe)r« 
\(N?Pg\Qg) T ,\ < c m ]T \g k \ \(N?Qg\e k )) r ,\ 

k 

< <?£ l<?fel |iVf Qg\r. < CX^ 2 £ \g k \{N^Qg\Qg) 1 ^ 



k k 

< C\- a ' 2 []T \g k f + m(N?Qg\Qg) r , 



< CX 



- a/2 [\Pg\^ + (N*Qg\Qg) 



r» 



where C > is a generic constant, which may differ from line to line. Hence for 
A > Ao, with Ao sufficiently large, we have 

(iVf g\g) r , = (N?Qg\Qg) r , + 2(N*Qg\Pg) Tm + (N*Pg\Pg) r , 



> \{N?Qg\Qg)T, + {N?Pg\Pg) T , - -fL\Pg\l m . 

This implies 



{B x g\g)T. = KTxg\g)v, + a(A*g\g)r, 

> c [^(iVf Qg\Qg)r, + A(JVf Pg\Pg)r,] 

+ c a{A*Qg\Qg) r ,-c\Pg\l t . 

Since is positive semidcfinitc and also A*Q has this property as im (Q) C 
-^2,o(r*), we only need to prove that X(N^ Pg\Pg)r t tends to infinity as A — > oo. 
To prove this, similarly as before we assume A > Ao and estimate 

|(iVf eile^r.)! < C|JVf < CX^ 2 {N» e^)*^. 

Choosing Ao sufficiently large this yields 

{N^Pg\Pg) L2(T ^ > c [min(7Vf e.leOL^p.) - S \Pg\l 2 (r,y 

Therefore it is sufficient to show 

lim X(N" e k \e k ) L2 , rr) = oo, fc = l,...,m. (3.17) 

A— >oo 

So suppose, on the contrary, that Xj(N^.g\g)i l2 (r t ) is bounded, for some g = e k 
and some sequence Aj — > oo. Then the corresponding solution i?j of (|3.13l) is such 
that := Aj^j is bounded in £2^) as 

Hence Uj has a weakly convergent subsequence, and we can assume without loss 
of generality that Vj — > Wqo weakly in Z,2(f2). Fix a test function ?/> £ T>(Q, \ T*). 
Then 

(K*Uj#)n = (d*Ai?j|V>)n = (#i|d*AV>)n = (vj|d„,A^) n /Aj ->■ 
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as j — y oo , hence Woo = in ^(fi). On the other hand we have 

0<HM = f g ds= [ -{dtdvtfjjds 
m Jr, Jr, 

= / d*A$j dx — Xj / K^flj dx — > / n*Vaodx, 
Jn Jn Jn 

hence Voo is nontrivial, a contradiction. This implies that f|3 . 1 T[) is valid. 

(d) Summarizing, we have shown that B\ is not positive semidefinite for small 
A > if cither T* is not connected or the stability condition does not hold, and 
B\ is always positive semidefinite for large A. Set 

Ao = sup{A > : is not positive semidefinite for each /i£ (0, A]}. 

Since B\ has compact resolvent, B\ has a negative eigenvalue for each A < Ao. 
This implies that is an eigenvalue of B\ , thereby proving that — L admits the 
positive eigenvalue Ao. 

Moreover, we have also shown that 

irj 

B h := lim XT\h + oA*h = c* 7 — — r— (/ — Qo)h + <rA*h. 

Therefore, B has the eigenvalue c*\T*\ /(nJl) L2 (n) ~cr(n— 1)/-R* with eigenfunc- 
tion e, and in case m > 1 it also possesses the eigenvalue — er(n — l)/R 2 with 
precisely (m — 1) linearly independent eigenfunctions of the form J2k a k^k with 
X^fc a k = 0- This implies that — L has exactly m positive eigenvalues if the stability 
condition does not hold, and m — 1 otherwise. 

(iv) (a) Suppose that (it, h) is an eigenfunction of L for the eigenvalue A = 0. 
Then dSHJ) yields 

2\^l /2 D\l + \dl /2 Vd\ 2 n = 0. (3.18) 

It follows from (|3.18l) and (13.6l) - (|3.7p that d is constant and D = on fl. Korn's 
inequality, in turn, implies Vu = on f2, and we then have u = by the no-slip 
condition on dfl. Moreover, the pressures are constant in the phases and we have 

[irj + <tAJi = 0, Ud-aA^h^Q on r*. 

We can now conclude from the relation 1*$ — crA^h = that the kernel of L is 
given by 

N(L) = span{(0, ^ , e), (0, 0, Y* ), . . . , (0, 0, F w fc ) : 1 < k < m}, (3.19) 

where the functions Yj = Yjfe k denote the spherical harmonics of degree one on 
r£, normalized by {Y*\ ^} fc )r fc = This shows that N(L) has dimension (mn+1), 
in accordance with the situation for the Stefan problem with surface tension [115] . 
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(b) It remains to show that A = is semi-simple if s ^ 0. We concentrate on 
the case where T* is connected, for simplicity. The disconnected case is treated in 
complete analogy. So suppose (u, &,h) € N(L 2 ). Hence L(u,$,h) £ N(L), i.e. 

n 

L(u, 0, h) = a (0, -a(n - l)/(l.I%),Y ) + £ a, (0, 0, Y), 

i=i 

where ao,a/ are appropriate coefficients and Yq = 1- Thus (it,$, h) solves the 
equations 



(3.20) 



and 





= 


in 


n\r, 


divu 


= 


in 


n\r* 


u 


= 


on 




M 


= 


on 


r. 


-[Tv] 


= — <t(A*/j)j/ 


on 





— d*Ai? = — aQK*(x(ri — in 17 \ 

= o on an 

[0] - on r» 

1*$ — aA*h = on T» 

t -(L/9*)u ■ v - {dMj = {h/e^^aiY on r., 



(3.21) 



We have to show ai = for all I. Integrating the equation for the temperature 
over fi we find 

a(n — !)(«;* | l)n '*|r*| 



a - 



a - 



(3.22) 



as u ■ v and the spherical harmonics Yi, 1 < I < n, all have mean zero on T*. 
Therefore, ao = 0, unless there is equality in the stability condition. If s ^ 0, and 
hence a = 0, it follows from (j3~2TJl) and (j3~2T|) 

= 2|^i /2 L>| 2 2 + a{A*h\u ■ v) Tt , 

n 

= 0»|di /2 V0| 2 2 - a(A*h\u ■ v + J2 a i Y i)r* = ft»|di /2 W|& - cr(A*% • v) Vt , 

;=i 

as A* is self-adjoint and AYi = for the spherical harmonics. Adding these 
equations gives 

2\n 1 J 2 D\l + Q*\d 1 J 2 \7d\l = 0. 

This implies D = 0, d constant, u ■ v — and m = 0, which in turn yields 
= Ym=i a iYi- Thus ai = for all I since the spherical harmonics Y. are linearly 
independent. Therefore, the eigenvalue A = is semi-simple. 

(v) Suppose for the moment that r» consists of a single sphere of radius R* — 
<j(n — 1)/ [?/>(#*)], centered at the origin of 1". Suppose S is a sphere that is 
sufficiently close to I\. Denote by (£i, ■ • ■ , £«) t ne coordinates of its center and let 
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Co be such that a(n — l)/\ip(9* + #*Co)] corresponds to its radius. We observe that 
the equation a(n — l)/[i/>(#* + #*Co)l = R has a unique solution Co for R close to 
i?*, as ^ by assumption. Then, by [HI Section 6], the sphere S can be 

parameterized over T* by the distance function 



p(0 = J2^- R * 

3=1 



\ 



n 

5 3 ' 

3=1 3=1 



(J2 CjYj) 2 + (a(n - + ^Co)]) 2 - 



Denoting by O a sufficiently small neighborhood of in R™ +1 , the mapping 

[C H- *(C) := (0,0.(1 + Co), p(C))] = O -> W*(nr +1 x W*-V*(T.) 
is C 1 (in fact C k if i/j is C k ), and the derivative at is given by 



*'(0)z = (0, 6*,-a(n - 1 WC^MW*)] 2 )*o + (0, 0, £ 

3 = 1 

Noting that 

<r(rc- 1)04^(0.)] = ^fl 2 
[V>(0*)] 2 a(n-l) 
we can conclude that near e» = (0, 0», T„) the set £ of equilibria is a C -manifold in 
W 2 (ft)™ x VF p 2 (fi) x VF p 4_1/p (r») of dimension (n + 1), and that T e ^ is isomorphic 
to the eigenspace N(L). 

It is now easy to see that this result remains valid for the case of m spheres of the 
same radius R„ . The dimension of £ is then given by (mn + 1), as mn parameters 
are needed to locate the respective centers, and one additional parameter is needed 
to track the common radius. □ 

Remarks 3.4. (a) One should observe that for the case s — 0, the eigenvalue A = 
ceases to be semi-simple and the dimension of the generalized eigenspace raises by 
one. This can be shown by similar arguments as in the proof of Theorem 2.1.(c) 
in [T7]. 

(b) For the Frechet derivative of the energy functional E(it, 0, T), see (|1.2p for the 
definition, we obtain 



ds . 



(E'(u, 0, r)|(«, 0, h)) = J " (u • v + e'{6)d) dx - j {aH T + \\\u\l + e(9)j)h 

At equilibrium (it, 9, T) = (0, 0*, r») this yields 

(E'(O,0„r»)|(«,0,fc)) = / n*$dx- [ {aHr,+{e4)hds. 

Here [ej := [6(61*)] = ty(0.)l - 9 m [if/(9 m )] = -(aH r , + h) where we used the 
equilibrium relation ["0(0*)] + crH-p, = and the definition of Z* in the last step. 
Preservation of energy then requires (E'(0, 0*,r*)|(w, d, h)) — 0, hence 

dx + I* J hds = 0. 
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In this case h = y^_ n QtjYi, where Yq — 1 and Yi denote the orthonormalized 
spherical harmonics of degree one. Hence h — a , and i3 = —a a(n — l)9*/l*B%, 
which implies 

~a(n — l)0*(«;*|l)n 



LRl 



0. 



Thus h = ao — unless we have equality in the stability condition (|3.1ip . Con- 
servation of energy kicks out one dimension of the eigenspace. 



4. Nonlinear Stability of Equilibria 

1. We now consider problem (|1.1[> in a neighborhood of a non-degenerate equi- 
librium e* = (0, 0*,r*) e £, with I* = 1(8*) ^ 0. Setting E = r* the transformed 
problem becomes 



dtu — /i,A« + V-7T 


= F u (u,ir,$,h) 


in 


n\s 


div u 


= F d (u,h) 


in 


fi \ s 


u = 0, d v ti 


= 


on 


an 


M=0, 


= 


on 


E 


-P e [m*(Vk + [Vm] t )K 


= G T (u,tf,h) 


on 


E 


-{T^-i^j+aAsh 


= G v (u,'9,h) 


on 


E 


Kxdt'd — g?*Ai9 


= F e (u,V,K) 


in 


q\e 




= G e (#,h) 


on 


E 


(i*/e*)(d t h-u-w)-ld*d u ti} 


= G h (u,d,h) 


on 


E 


u(0) = u , 0(0) = 0o, Mo) 


= ft-o- 







(4.1) 



The nonlinearities on the right hand side of (14.11) are, up to some straightforward 
modifications, defined in Section 7}. It follows that the nonlinearities are of 
class C 1 from E( J) to F(J), and they satisfy 

^(0) = G fc (0)=F;(0) = G^(0) = 0, je{ u ,d,n, fce{T,M,M- 

In order to shorten notation, we will occasionally write (|4.1[) in short form 

hz = N(z), z(Q) = z . 

The state manifold locally near the equilibrium = (0, T*) reads as 

SM := Uu,#,h) G (% 2 ~ 2 / p (fi\E)nC(fi)) n+1 x iy p 4 " 3 /P(E), 

div u = F d (u, h) in Q \ E, u = 9„0 = on 90, (4.2) 
-P s [MVm+ [Vti]> s = G T (u,d,h) on E, 
i,0 - <7A E /i = G e (0, /i), [d»0 v 7?] + G h (u, i?, h) e V^- 6 / p (E)}. 
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Note that due to the compatibility conditions this is a nonlinear manifold. We 
shall parameterize this manifold over its tangent space 

Z :=Uu,$,h) € (W^ 2/p (!l\S)nC(li))" +1 x iy p 4 - 3/p (S), 

div u = in i7 \ E, u — d^tf — on dfl, 

-P s [^(Vu+[Vu] T ]^ s =0, lJ-<rAvh = 0, [d,0„0] € W^ 6 / p (E)}. 
We mention that the norm in Z is given by 

= N^-a/f^E) + l^l Wp 2 - 2 /P(f2\ E ) + Nw* _3/:P (E) + I M*^] I lV p 2 " 6/p (S) ' 

2. In order to parameterize the state manifold SAA over Z near the given equi- 
librium (0, (9», E) we consider the linear elliptic problem 

LUU — /i*Alt + V7T = 

div u = / d 
u = 0„i? = 

[u] = M = o 

-k(V« + [Vti]> s + [tt]i/e + o-(A s /i>£ - 3„ 

— d* A?9 = 

— aA^h = go 
(U/0*)(ujh - u ■ v?) - \d*d v -d] = g h 

for given data (fd, 9u, 9h, 9e)- For this problem we have the following result. 

Proposition 4.1. Suppose p > 3, i* =^ and w > is sufficiently large. Then 
problem (|4.3p admits a unique solution (u, tt, i9, /i) wf/j regularity 

(M,fc)e (% 2 - 2/p (n\s)nc(!i))'' +1 xT^ p 4 - 3 / p (s), 7 re% 1 - 2 / p (r!\s), 

if and only if the data (fd, g u , gh, ge) satisfy 

fd G W^ 2 ^ \ E) n ( 5u , 5 ,0 G ^-^(E^ 1 , 5e G ^ 2 - 3 / p (E). 

T/ie solution map [{fd,gu,gh,ge) | — > (u, 7r, 1?, /i)] is continuous in the corresponding 
spaces. 

Proof. This purely elliptic problem can be solved in the same way as the corre- 
sponding linear parabolic problem. □ 

Theorem 4.2. There exists a neighborhood U of in Z and a map 
dl> G C X {U, (W^- 2/p (n \ E) n x W p 4 - 3 /f(E)) with 0(0) = 0'(O) = 0, 

such that [z 4 z + <K^0] ■ ^ ~ ^ 5.M provides a parameterization of the state 
manifold SA4 near the equilibrium (0, 0*, E). 



111 
in 
on 
on 
on 
in 
on 
on 



n\ e 
n\ e 

E 
E 

E 
E 



(4.3) 
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Proof. Fix any large u) > 0. Given z = (u, $,h) € Z sufficiently small, and setting 
(u, $, h) = (u, h) + (u, h), we solve the nonlinear elliptic problem 



UJU — [1,Am + V7T 


= 


in 


n\s 


divu 


= F d (u,h) 


in 


n\s 


u — d v d 


= 


on 


an 


m = m 


= 


on 


s 


-P s [ M *(Vw+[Vz2] T ]^ 


= G T (u,$,h) 


on 


E 


-([M*(VG + [VS] T K|^) + [tt] + aAvh 


= 


on 


S 




= 


in 


n\s 


h$ — aA^h 


= G g {0,h) 


on 


s 




= G h (u,ti,h) 


on 


s 



by means of the implicit function theorem, employing Proposition 14. f I Then with 
z = (u,$,h) and z = z + z we obtain z — 4>{z), with a C^-function <fi such that 
0(0) = 0'(O) = 0. Then z = z + <p(z) E SA4, hence SM. is locally parameterized 
over Z . 

To prove surjectivity of this map, for given (u, i?, h) G SA4, solve problem (|4.4I) . 
where the functions (Fd(u, h), G T (u, i?, h),Gg($, h), Gh(u, hj) are now given. By 
Proposition 14.11 the resulting linear problem has a unique solution z = (u, 
Let z = z — I. Then we see that z = </>(z), hence the map [z M> z + 0(5)] is also 
surjective near 0. □ 



3. Next we derive a similar decomposition for the solutions of problem (|4.1I) . 
Let z = (z o ,0(zo)) G S.M be given, and let z € E(J) be the solution of (I4.1[) 
with initial value zo. Then we would like to devise a decomposition of z such 
that z(t) = z(t) + z(t) with z(t) G Z for t > 0. As before, we use the notation 
z = (u, /i), and z = (u, i9, h). In order to accomplish this we consider the coupled 
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systems of equations 



LOU + 8tU — fJ,*Au + V7T 


= F u (u,ir,{),h) 


in 


0\E 


div u 


= F d (u,h) 


in 


n\s 


u — d v -d 


= 


on 




[uj = p] 


= 


on 


s 


-P s [/^(Vu+[Vu] T )]^ 


= G T (u,$,h) 


on 


s 


-([M*(Vu+ [Vm] t 1^ e |^e) + [#] + (T A s h 


= G„{u,$,h) 


on 


s 


K*Lud + n^dtfi — d*A?9 


= F e {u 1 d 1 h) 


in 


n\s 


l*"d — aAsh 


= G e {$,h) 


on 


s, 


{h/6*)uh + (l*/6*)(d t h -u-v?)- [d*«9„tf] 


= G h (u,#,h) 


on 


s 


2(0) 


= <t>{z ) 







(4.5) 



and 



<9 t U — fl^Au + V-7T 


= UJU 


in 


n\s 


div u 


= 


in 


n\s 


u = du'd 


= 


on 


50 


[«] = M 


= 


on 


E 


-ft[/4(V5+[VS]> s 


= 


on 


S 


([/i*(Vu + [Vu] T ]i/ s |z/ s ) + [ttJ + crA E ft 


= 


on 


S 


K*dt& - d*Atf 




in 


n\s 


— a A-^h 


= 


on 


s, 


{h/6*)(d t h-u-vx)-ld*dj} 


= (l*/6*)ujh 


on 


E 


m 


= z . 







(4.6) 



Equations (|4. 5|) (|4.6p can be rewritten in the more condensed form 

L u z = N(z + z), z(O) = 0(z o ) 
z + Lz = wz, 5(0) = Sq. 



(4.7) 



where we use the abbreviation to denote the linear operator on the left hand 
side of (|4.5p . and N to denote the nonlinearities on the right hand side of (|4.5[) . 
respectively. 

Remark 4.3. As divu is in general nonzero, z does not belong to the base 
space Xq- However, this defect can be easily overcome, replacing u in (|4.6p by 
its Hclmholtz-projection in Q. This only changes the pressure tt by a jump-free 
part, but the velocity u, and hence also z, are unchanged. 
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4. For the purpose of proving the stability result it turns out to be more convenient 
to modify the decomposition of z derived in the previous step in the following way. 
Suppose = Zoo + 4>(zoo) € SM. is an equilibrium of (I4.1[) which is close to the 
fixed equilibrium z* = (0, 9*,T*). Then we decompose the solution z of (|4.1[) as 
z(t) = Zoo + z(t) + z(t), where as above z(t) £ Z. Clearly, Lzoo = N(z OQ ), and we 
are lead to consider the following coupled system for the pair (z, z) 

h u z = N(z QO + z + z)- N(zoo), z(0) = 4>(z Q ) - (f>(Zoo), 

(4.8) 

z + Lz = ujz, 5(0) = zq — Zqo. 



The abstract problem (|4.8p can be treated in the same way as in the proof of 
Theorem 5.2 in Priiss, Simonett and Zacher [19]. This implies the following result. 

Theorem 4.4. Let p > n + 2, a > and £* ^ 0, and suppose tpi £ C 3 (0, oo), 
^i,di £ C 2 (0, oo) are such that 

Ki(s) = -stp'/(s) > 0, Hi(s) > 0, di(s) > 0, s£ (0, oo), j = 1, 2. 

Le< the function ip be as in (|3.2p . T/ien m i/ie topology of the state manifold SA4 
we have: 

(a) (O,0*,r„) € £ is stable ifT* is connected and ip'{6*) < 0. 

Any solution starting in a neighborhood of such a stable equilibrium con- 
verges to another stable equilibrium exponentially fast. 

(b) (0,0*, r*) G £ is unstable i/T* is disconnected or tp'(9*) > 0. 

Any solution starting and staying in a neighborhood of such an unstable 
equilibrium converges to another unstable equilibrium exponentially fast. 

5. Global Existence and Convergence 

We have seen in [16] that the negative total entropy, see (|1.3p . is a strict Lya- 
punov functional. Therefore the limit sets of solutions in the state manifold SM 
are contained in the manifold £ C SM. of equilibria. 

There are several obstructions against global existence: 

• Regularity: the norms of either u(t), 9(t), T(t), or [d(#(<))<9„#(i)] may 
become unbounded; 

• Well-posedness : the condition 1(9) ^ may be violated; or the tempera- 
ture may become 0; 

• Geometry: the topology of the interface may change; 
or the interface may touch the boundary of Q; 

or a part of the interface may contract to a point. 
Recall that the compatibility conditions 

divu(t) = in Q \ ]?(£), u(t) = d u 6(t) = on dfl, 
[u(t)] = p] = P r M9(t))D(t)] = 0, + <7flr(t) = on T(t), 

are preserved by the semiflow. 
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Let (u, 9, r) be a solution in the state manifold SAi with maximal interval of 
existence [0, £*). By the uniform ball condition we mean the existence of a radius 
7'o > such that for each t £ [0,t*) ; at each point x £ T(t) there exists centers 
Xi £ ili(t) such that B ro (xi) C fli and T(t) D B ro (xi) = {x}, i = 1,2. Note that 
this condition bounds the curvature of T(t), prevents parts of it to shrink to points, 
to touch the outer boundary <9f2, and to undergo topological changes. 

With this property, combining the local semiflow for (jl.ljl with the Lyapunov 
functional and compactness we obtain the following result. 

Theorem 5.1. Letp > n+2, a > 0, and suppose ipi £ C 3 (0,oo), fJ>i,di £ C 2 (0,oo) 
such that 

Ki(s) = -sip"(s) > 0, fii(s) > 0, di(s) > 0, s £ (0, oo), i = 1,2. 

Suppose that (u, 0, T) is a solution of in the state manifold SAi on its maxi- 

mal time interval [0,i*). Assume there is constant M > such that the following 
conditions hold on [0,t*): 

(i) \u(t)\ [w i- VP]n , \9(t)\ w i- Vp , |r(t)| w 4- 3/P , \ld(9(t))d v 6(t)]\ w 2- e/P < M; 

(ii) \l(9(t))\,9(t)>l/M; 

(hi) r(t) satisfies the uniform ball condition. 
Then t* = oo, i.e. the solution exists globally, and its limit set w(it,#,r) £ £ is 
non-empty. If further (O^ocToo) £ 9,T) with Too connected and V5'(#oo) < 

0, then the solution converges in SA4 to this equilibrium. 

Conversely, if (u(t), 9(t), T(t)) is a global solution in SA4 which converges to 
an equilibrium (0, 0*, T„) £ £ in SAi as t — > oo, and 1(9*) ^ 0, then (i)-(iii) hold. 

Proof. Under the assumptions (i)-(iii) it is shown in the proof of |161 Theorem 
8.2] that = oo and that the orbit (u, 9, T)(]R + ) C SAi is relatively compact. 
The negative total entropy is a strict Lyapunov functional, hence the limit set 
uj(u, 9,T) C SAi of a solution is contained in the set £ of equilibria. By com- 
pactness, w+(u,#,r) C SAi is non-empty, hence the solution comes close to £, 
and stays there. Then we may apply the convergence result Theorem 14.41 The 
converse follows by a compactness argument. □ 

Remarks 

(i) We believe that in Theorem 15. II the assumption that is connected can be 
dropped and ^'(^oo) < can be replaced by </?'(#oo) 7^ 0. However, a proof of this 
requires much more technical efforts, we refrain from these, here. 

(ii) We cannot show that the temperature stays positive if it is initially since we 
did not make any assumptions on the behavior of the functions /ij, dj, ipj near 0. 
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